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1.  INTRODUCTION 


In  2003,  the  Defense  Seienee  Board  observed:  “The  ...  problem  is  that  instruments  that  ean 
deteet  the  buried  UXOs  also  deteet  numerous  serap  metal  objeets  and  other  artifaets,  whieh  leads 
to  an  enormous  amount  of  expensive  digging.  Typieally  100  holes  may  be  dug  before  a  real 
UXO  is  unearthed!  The  Task  Foree  assessment  is  that  mueh  of  this  wasteful  digging  ean  be 
eliminated  by  the  use  of  more  advaneed  teehnology  instruments  that  exploit  modern  digital 
proeessing  and  advaneed  multi-mode  sensors  to  aehieve  an  improved  level  of  diserimination  of 
serap  from  UXOs.”  In  keeping  with  these  remarks  and  with  prior  funding  (UX-1225,  MM-0437, 
and  MM-0838),  the  LBNL  group  has  suoeessfully  designed  and  built  the  eart-mounted  Berkeley 
UXO  Diseriminator  (BUD)  and  demonstrated  its  performanee  at  various  test  sites  (ef. 
Gasperikova  et  ah,  2007,  2008,  and  2009). 

Beeause  hand-held  systems  have  the  advantage  of  being  lightweight,  eompaet,  portable,  and 
deployable  under  most  site  eonditions,  they  are  partieularly  useful  in  areas  of  dense  vegetation  or 
ehallenging  terrain.  In  heavily  wooded  areas  or  areas  with  steep  or  uneven  terrain,  hand-held 
sensors  may  be  the  only  suitable  deviee  for  UXO  deteetion  and  diserimination  beeause  it  ean  be 
earned  through  spaees  that  the  operator  eould  walk  through  or  at  least  approaeh.  Furthermore,  it 
is  desirable  to  find  and  eharaeterize  a  metallie  objeet  without  the  need  to  aeeurately  loeate  the 
sensors  at  multiple  positions  around  the  target.  The  ideal  system  would  thus  loeate  and 
eharaeterize  the  target  from  a  single  position  of  the  sensor  and  indieate  to  the  operator  where  to 
flag  the  target  for  subsequent  study. 


4 


Based  on  these  eonsiderations,  we  designed  and  built  a  sensor  paekage  in  a  shape  of  a  14-in 
(0.35  m)  eube.  This  hand-held  prototype  incorporates  the  key  features  of  the  cart-mounted 
system  -  (a)  three  orthogonal  transmitters  and  ten  pairs  of  receivers,  and  (b)  difference  or 
gradient  measurements  that  significantly  reduce  the  ambient  and  motion  noise,  and  greatly 
enhance  the  sensitivity  to  the  gradients  of  the  target.  The  system  characterizes  the  target  from  a 
single  position.  Results  from  a  local  test  site  were  in  a  good  agreement  with  theoretical 
performance  calculations  of  such  a  device.  This  survey  was  designed  to  demonstrate 
performance  of  the  system  under  realistic  survey  conditions  at  the  Aberdeen  Proving  Ground 
(APG)  in  Aberdeen,  MD.  The  survey  was  preformed  in  an  area  with  known  items  (“Calibration 
Grid”),  and  in  a  seeded  blind  test  area  (the  “Blind  Test  Grid").  The  ground  truth  for  the  surveys 
conducted  on  the  Blind  Test  Grid  was  withheld  from  the  testers.  Only  the  graded  test  scores 
based  on  target  detection  and  target  classification  were  provided.  For  more  information,  see 
http;//aec. army.mil/usaec/technology/uxoO  lc03.html. 
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2,  TECHNOLOGY  DESCRIPTION 


The  SERDP  and  ESTCP  have  supported  EBNL  in  development  of  a  deviee  that  not  only  deteets 
the  objeet  itself  but  also  quantitatively  determines  its  size,  shape,  and  orientation.  In  a  hand-held 
system  the  sensor  paekage  must  be  easily  maneuvered  over  rough  terrain,  through  brush,  ete.  A 
useful  eriterion  might  be  that  it  eould  be  earned  through  spaees  that  the  operator  eould  traverse  or 
at  least  approaeh.  To  avoid  the  relianee  on  aeeurate  multiple  positioning  of  any  system  it  has  been 
shown  that  (1)  the  objeet  must  be  illuminated  with  three  different  polarizabilities  of  the  transmitted 
field  and  that  (2)  to  determine  its  loeation  from  a  single  position  of  the  transmitter-reeeiver  system 
multiple  reeeivers  must  be  used.  To  aoeommodate  the  first  requirement  a  hand-held  design 
ineorporates  three  orthogonal  transmitter  loops  mueh  like  the  eart-mounted  system.  To 
aoeommodate  the  seoond  requirement  the  hand-held  UXO  disoriminator  uses  ten  pairs  of  reeeivers. 
The  same  disorimination  abilities  afforded  by  the  eart-mounted  system  are  available  in  the  hand¬ 
held  unit  albeit  with  slightly  reduced  depth  of  detection.  The  hand-held  UXO  discriminator  is  able 
to  discriminate  small  (20  mm)  objects  at  a  depth  of  0.45  m  and  large  (105  and  155  mm)  objects  at  a 
depth  of  0.85  m  and  to  detect  them  down  to  1.15  m.  The  hand-held  prototype  is  shown  in  Eigure  1. 
The  transmitter-receiver  configuration  is  shown  in  Eigure  la  while  the  assembled  transmitter- 
receiver  cube  is  shown  in  Eigure  lb. 
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Figure  1,  LBNL  Hand-held  UXO  Discriminator  -  (a)  transmitter-receiver  array  configuration, 
(b)  assembled  prototype. 


The  hand-held  UXO  discriminator  employs  three  orthogonal  transmitters  and  ten  pairs  of 
differenced  receivers.  Each  vertical  face  of  the  cube  has  three  induction  coils,  while  the  two 
horizontal  faces  have  four  induction  coils,  each  sensitive  to  the  magnetic  field  component  normal 
to  the  face.  Receivers  on  opposite  faces  of  the  cube  are  paired  along  the  symmetry  lines  through 
the  center  of  the  system  and  each  pair  sees  identical  fields  during  the  on-time  of  current  pulses  in 
the  transmitter  coils.  They  are  wired  in  opposition  to  produce  zero  output  during  the  on-time  of 
the  pulses  in  three  orthogonal  transmitters.  This  configuration  dramatically  reduces  noise  in 
measurements  by  canceling  background  electromagnetic  fields  (these  fields  are  uniform  over  the 
scale  of  the  receiver  array  and  are  consequently  nulled  by  the  differencing  operation),  and  by 
canceling  noise  contributed  by  the  tilt  of  the  receivers  in  the  Earth’s  magnetic  field,  thus  greatly 
enhances  receivers’  sensitivity  to  gradients  of  the  target  response. 

The  transmitter  coils  are  powered  separately  from  the  data  acquisition  board.  Pulsers  provide 
resonant  circuit  switching  to  create  bi-polar  half-sine  pulses  of  300  ps  width.  The  current  has  a 
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~40  A  peak  and  a  resonant  eireuit  voltage  of  -400  Volts.  The  operational  overall  half-sine  duty 
eyele  is  -12%.  Transients  are  digitized  with  a  sampling  interval  of  4  ps.  The  sensors  are 
eritieally  damped  5-meh  eoils  with  a  self-resonant  frequeney  of  75  kHz.  Data  aequisition  is 
performed  on  a  single  board.  This  board  has  12  high-speed  ADC  ehannels  for  the  output.  Ten  of 
these  ehannels  are  used  for  the  signal  from  reeeiver  eoils,  and  the  remaining  two  ehannels 
provide  information  about  the  system  (i.e.  tilt  information,  transmitter  eurrent). 

An  important  feature  of  the  hand-held  UXO  diseriminator  is  an  inversion  algorithm,  whieh  is  used 
to  determine  target  properties  from  measurements  with  a  given  transmitter-reeeiver  eonfiguration. 
At  any  given  time  the  response  is  inverted  to  yield  the  loeation  (x,  y,  z)  of  the  target,  its  attitude 
and  its  prineipal  polarizabilities,  whieh  yield  an  apparent  aspeet  ratio.  Signal-to-noise  (S/N) 
estimates  (or  measurements)  are  interpreted  in  this  inversion  to  yield  error  estimates  on  the 
loeation,  attitude  and  polarizabilities. 

It  has  been  demonstrated  that  a  satisfaetory  elassifieation  seheme  is  one  that  determines  the 
prineipal  dipole  polarizabilities  of  a  target  -  a  near  intaet  UXO  displays  a  single  major 
polarizability  eoineident  with  the  long  axis  of  the  objeet  and  two  equal  transverse  polarizabilities. 
The  indueed  moment  of  a  target  depends  on  the  strength  of  the  transmitted  indueing  field.  The 
moment  normalized  by  the  indueing  field  is  the  polarizability.  This  deseription  of  the  inherent 
polarizabilities  of  a  target  is  a  key  in  diseriminating  UXO  from  irregular  serap  metal.  Figures  2 
and  3  illustrate  a  diserimination  eapability  of  the  system  for  UXO  objeet  and  a  serap  metal, 
respeetively.  While  UXO  objeets  have  a  single  major  polarizability  eoineident  with  the  long  axis 
of  the  objeet  and  two  equal  transverse  polarizabilities  (Figure  2),  serap  metal  exhibits  three 
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distinct  principal  polarizabilities  (Figure  3).  There  are  clear  distinetions  between  symmetrie 
intact  UXO  and  irregular  scrap  metal.  Moreover,  UXO  have  unique  polarizability  signatures,  and 
thus  distinetions  ean  be  made  among  various  UXOs. 


Object  orientation  and  equivalent  dipole  polarizability  estimates  of  large  and  shallow  UXO  or 
scrap  are  more  problematie  as  they  are  affected  by  higher  order  (non-dipole)  terms  induced  in 
objeets  due  to  souree  field  gradients  along  the  length  of  the  objects.  In  the  case  of  the  large  and 
shallow  objects,  the  hand-held  system  can  be  easily  raised  an  appropriate  distance  above  the 
object  such  that  the  dipole  model  approximation  for  polarizability  estimates  is  appropriate. 
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Figure  2,  Inversion  results  for  the  principal  polarizabilities  of  81  mm  mortar. 
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Figure  3,  Inversion  results  for  the  prineipal  polarizabilities  of  19x8  em  serap  metal. 
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3.  DEMONSTRATION  SURVEY 


The  demonstration  survey  was  performed  between  June  21,  2010  and  July  1,  2010.  The  hand¬ 
held  UXO  diseriminator  was  operated  in  a  cued  mode.  The  system  was  brought  to  marked 
locations  and  run  in  the  characterization/discrimination  mode.  The  three  discriminating 
polarizability  responses  along  with  the  object  depth  and  horizontal  location  with  respect  to  the 
center  of  the  bottom  plane  of  the  cube  were  recorded  and  visually  presented  on  the  computer 
screen.  Additional  values  recorded  at  each  location  were  the  S/N  ratio,  and  a  polarizability  index, 
which  is  an  average  value  of  the  product  of  time  (in  seconds)  and  polarizability  rate  (in  m  /s) 
over  the  46  sample  times  logarithmically  spaced  from  80  to  1460  p,s.  The  survey  took  two 
weeks.  The  schedule  accommodated  partial  days  when  the  base  was  closed.  We  had  a  one-day 
delay  because  of  a  problem  with  the  system.  An  unusual  and  unexpected  heat  wave  made  the 
first  half  of  the  survey  extra  challenging.  The  system  performed  well,  and  we  collected  a  high 
quality  data.  Figure  4  illustrates  the  deployment  of  the  hand-held  UXO  discriminator  in  the  field. 


Figure  4,  LBNL  hand-held  UXO  discriminator  at  APG. 
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The  Calibration  Grid  consists  of  66  cells  4  m  on  a  side  arranged  in  11x6  layout  (A-K  and  1-6). 
The  APG  personnel  marked  all  66  cell  centers  locations  with  plastic  pin  flags.  These  cells 
contain  14  different  munitions  types  (from  20  mm  projectile  to  155  mm  projectile),  and  eight 
representative  clutter  items  buried  at  various  orientations  and  depths.  The  shallowest  object  is  a 
horizontal  40  mm  projectile  at  0.0  m  depth,  and  the  deepest  object  is  a  horizontal  155  mm 
projectile  at  1.08  m  depth.  We  took  measurements  along  lines  2  m  apart,  with  a  2  m  distance 
between  the  measurements  on  each  line.  In  addition  to  soundings  at  centers  of  grid  cells,  we 
collected  data  at  0.15  m  before  and  0.15  m  after  cell  centers  along  survey  lines.  When 
polarizability  inversions  didn’t  agree  in  appearance,  provided  that  the  central  sounding  had  S/N 
larger  than  four,  eight  additional  soundings  were  acquired,  at  0.3  m  and  0.5  m  offset  on  the  arms 
of  a  cross,  centered  on  the  cell  center. 

The  Blind  Test  Grid  is  40x40  m  area.  The  individual  cells  are  1x1  m  with  an  empty  1x1  m  cells 
in  between.  There  are  400  opportunities/flags  arranged  in  20x20  layout  (A-T  and  1-20).  Again, 
all  400  points  were  marked  with  plastic  pin  flags.  The  Blind  Test  Grid  is  made  up  of  six  different 
types  of  munitions  and  clutter  found  in  the  Calibration  Grid,  and  in  the  direct  and  indirect  fire 
sub-areas  of  the  Open  Field.  We  took  measurements  along  the  lines  1  m  apart,  with  aim 
distance  between  the  measurements  on  each  line.  In  addition  to  soundings  at  the  test  grid  cell 
centers,  we  collected  data  at  0.15  m  before  and  0.15  m  after  cell  centers  along  the  survey  lines. 
When  polarizability  inversions  didn’t  agree  in  appearance,  provided  that  the  central  sounding 
had  S/N  larger  than  four,  eight  additional  soundings  were  acquired,  at  0.3  m  and  0.5  m  offset  on 
the  arms  of  a  cross,  centered  on  the  cell  center. 
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4,  DATA  ANALYSIS  AND  INTERPRETATION 


Ten  channels  of  field  data  were  recorded  at  a  rate  of  250  k-samples/second  for  each  of  three 
transmitters.  Field  data  were  stacked  together  in  a  field  programmable  gate  array  (FPGA)  and 
transferred  to  a  field  computer  (laptop)  forming  a  primitive  stack.  An  even  number  of  primitive 
stacks  was  averaged  together  to  form  stacked  data  for  further  processing.  The  peak  transmitter 
current  was  estimated  from  the  stacked  transmitter  current  channel  record,  and  the  data  were 
normalized  by  that  value.  Nominal  transmitter  shut-off  time  was  estimated,  and  induction 
responses  were  computed  at  46  logarithmically  spaced  times  between  80  and  1460  ps,  averaged 
in  half-sine  windows  with  widths  10%  the  center  time  after  transmitter  pulse  shut-off  Responses 
were  differenced  with  background  responses  collected  over  an  area  determined  to  be  free  of 
metallic  objects.  The  resulting  30  channels  of  normalized  responses  were  then  inverted  for  a 
candidate  object  position  and  principal  polarizabilities  as  a  function  of  time  after  transmitter 
shut-off. 


4.1  Training  Data 

The  empirical  densities  for  UXO  and  scrap  classes  are  derived  from  the  training  data.  The  UXO 
class  training  data  included  polarizability  responses  collected  over  UXO  at  the  local  test  facility 
as  well  as  at  the  APG  Calibration  Grid.  Both  sets  contain  responses  from  14  different  munition 
types  (from  20  mm  projectile  to  155  mm  projectile)  buried  at  various  orientations  and  depths. 
These  responses  are  considered  reliable  at  observation  times  from  80  ps  to  1460  ps  after  the 
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extinction  of  the  primary  field.  These  were  supplemented  with  an  approximately  equal  number  of 
additional  UXO  polarizability  responses  previously  collected  using  the  cart-mounted  system 
(BUD).  Adding  them  to  the  training  data  allowed  for  better  estimates  of  the  densities  of 
polarizability  responses.  For  the  scrap  response  class,  the  responses  from  the  APG  Calibration 
Grid  were  supplemented  with  BUD  scrap  responses  from  previous  surveys  at  Camp  Sibert,  AL, 
and  San  Luis  Obispo,  CA.  The  BUD  responses  are  considered  reliable  from  140  ps  to  1400  ps. 
To  extend  them  to  the  window  of  the  hand-held  prototype  responses,  least  squares  predictions  of 
responses  and  suitable  noise  was  based  on  the  responses  from  the  original  BUD  time  window, 
with  prediction  coefficients  based  on  the  hand-held  responses. 


4.2  Calibration  Grid 


The  Calibration  Grid  is  44x24  m  area  with  cells  in  11x6  layout.  All  66  centers  were  marked  with 
plastic  flags,  since  at  the  moment  our  system  doesn’t  have  GPS  navigation  capabilities.  We  took 
measurements  along  lines  2  m  apart,  with  a  2  m  distance  between  the  measurements  on  each  line. 
In  addition  to  the  soundings  at  the  grid  cell  centers,  we  collected  data  at  0.15  m  before  and  0.15  m 
after  the  cell  center  along  the  survey  lines.  The  soundings  were  differenced  with  the  background 
reference  soundings  taken  within  the  previous  30-40  minutes  at  a  nearby  site  determined  by  the 
field  operator  to  be  free  from  metallic  objects.  An  estimate  of  S/N  was  made  based  on  signal 
levels  relative  to  median  drift  amounts  recorded  during  the  tests  at  our  field  test  facility  in 
California  (Richmond  Field  Station,  RFS).  The  estimate  is 
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where  p(r)  is  the  negative  of  a  log  likelihood  function  estimated  from  median  normalized  system 
drifts  between  reference  soundings  taken  before  and  after  UXO  soundings,  p(r)''  is  its  inverse, 
Sijk  is  the  signal  at  the  i'th  time  in  the  j'th  receiver  pair  in  the  response  to  the  k'th  transmitter, 
medianijk  is  the  median  absolute  value  of  the  drift  for  the  ijk'th  time  receiver  pair  and  transmitter 
combination  in  the  stability  measurements,  N  is  the  product  of  the  number  of  receiver  pairs  (10), 
the  number  of  transmitters  (3),  and  the  number  of  response  time  averages  considered  (46),  which 
is  used  also  for  the  summation. 

The  Calibration  Grid  S/N  map  is  shown  in  Figure  5.  The  recovered  responses  (66)  of  the  objects 
at  the  center  of  each  cell  were  used  as  the  training  data  for  the  Blind  Test  Grid  objects 
discrimination  and  identification.  The  polarizability  plots  of  these  objects  are  shown  in  Appendix 
1 .  The  title  of  the  plot  consists  of  a  prefix  CAL  indicating  that  data  are  from  the  Calibration  Grid, 
followed  by  the  x  cell  location  (x  goes  from  0  to  44  m,  from  left  to  right,  which  is  a  west-east 
direction  in  Figure  5),  followed  by  the  y  cell  location  (y  goes  from  0  to  24  m,  from  top  to  bottom, 
which  is  a  north-south  direction  in  Figure  5).  The  rest  of  the  file  name  can  be  ignored  -  Run  and 
four-digit  number  is  a  counter  that  guarantees  uniqueness  of  the  file  name,  and  various 
extensions  between  periods  indicate  a  sequence  of  the  files  used  in  the  inversion.  For  example, 
CAL-x34-y22-Run-0346.bin-diff.avg.inv  is  a  response  of  the  cell  with  x=34  m  and  y=22  m, 
which  is  the  cell  16.  We  took  additional  sub-grid  measurements  around  the  cell  center  when 
recovered  polarizabilities  significantly  departed  from  a  nominal  response  of  the  object  in  that  cell 
or  S/N  was  too  low.  Based  on  our  detection  threshold  of  S/N  less  than  four,  cell  C3  would  be 
considered  empty,  while  in  reality  it  contained  a  small  piece  of  scrap.  Cells  with  S/N  between 
four  and  six,  which  were  considered  as  “can’t  analyze”  in  our  analysis,  were  H4  (60  mm  mortar 
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at  0.6  m  depth)  and  K4  (105  mm  projectile  at  1.05  m  depth).  Cell  A6  contained  the  deepest 
object  -  155  mm  projectile  at  1.08  m  depth.  Based  on  our  previous  measurements  this  would 
produce  a  low  S/N.  However,  because  of  a  surface  nearby  metallic  object,  this  cell  S/N  was  high, 
although  we  were  not  able  to  identify  this  deep  object  correctly. 

If  one  compares  the  object  depths  estimated  by  the  inversion  algorithm  and  the  true  object 
depths,  56%  of  the  object  depths  were  within  ±0.05  m,  33%  were  between  0.05  and  0.15  m,  and 
11%  had  the  depth  differences  more  than  0.15  m.  The  cells  with  the  largest  difference  in  the 
object  depth  were  E5,  G3,  J2,  and  J3.  Cells  A6,  H4,  and  K4  that  we  could  not  analyze  also  fall 
into  this  category. 

Table  1  shows  the  Calibration  Grid  discrimination  results.  The  first  column  is  the  cell  name;  the 
second  column  is  the  calculated  probability  of  being  UXO  (p(UXO))  (anything  with 
p(UXO)>0.0005  was  considered  as  UXO,  9999  indicates  data  in  “can’t  analyze”  category);  the 
third  column  is  the  estimated  object  depth;  the  forth  column  is  the  object  identification  by  our 
algorithm;  the  fifth  column  is  the  ground  truth;  the  sixth  column  is  the  true  object  depth;  the 
seventh  column  is  the  file  used  for  the  analysis  (this  is  the  same  as  the  plot  title  in  the  Appendix 
1). 
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Figure  5,  Signal-to-noise  ratio  map  of  the  APG  Calibration  Grid.  Colors  are  plotted  on  a 
logarithmic  scale  for  better  visibility.  The  detection  threshold  of  S/N=4  is  equal  to  0.6  on 
the  logarithmic  scale. 


4.3  Blind  Test  Grid 


The  soundings  were  collected  at  the  test  grid  cell  centers  (400  opportunities),  and  at  0.15  m 
before  and  0.15  m  after  the  cell  center  along  the  survey  lines,  with  the  system  oriented  in  a  single 
direction.  The  soundings  were  differenced  with  the  background  reference  soundings  taken  within 
the  previous  30-40  minutes  at  a  nearby  site  determined  by  the  field  operator  to  be  free  from 
metallic  objects.  The  estimate  of  S/N  was  calculated  as  described  above.  The  cells  with  S/N 
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values  greater  than  four  were  eonsidered  oecupied.  When  S/N  values  were  over  400  the  system 
was  either  raised  0.2  m  or  0.4  m,  or  moved  away  by  0.25  and  0.5  m  to  reduee  the  S/N  below  400 
and  the  measurements  were  repeated. 

When  polarizability  inversions  from  soundings  taken  at  a  eommon  height  agreed  in  appearanee, 
then  the  sounding  nearest  laterally  to  the  interpreted  objeet  was  used  for  diserimination. 
Otherwise,  provided  that  the  central  sounding  had  S/N  larger  than  four,  eight  additional 
soundings  were  acquired,  at  0.3  m  and  0.5  m  offset  on  the  arms  of  a  cross,  centered  on  the  cell 
center.  When  polarizability  inversions  from  the  soundings  agreed  in  appearance  then  the 
sounding  nearest  to  the  interpreted  object  was  used  for  discrimination.  When  a  subset  of 
inversions  agreed  in  appearance  then  the  subset  member  closest  to  the  interpreted  object  was 
used  for  discrimination.  When  polarizability  inversions  from  all  soundings  in  a  cell  were  of 
disparate  appearance,  multi-site  inversion  was  performed  on  the  soundings,  and  if  it  fit  the  data 
to  better  than  50%  at  the  site  closest  to  the  interpreted  object,  the  resultant  dipole  polarizabilities 
were  used  for  discrimination,  otherwise  data  from  this  cell  were  considered  as  “can’t  analyze”. 
Similarly,  any  sounding  with  an  estimated  S/N  below  six,  which  our  discrimination  procedure 
estimated  to  be  scrap,  was  considered  as  “can’t  analyze”.  The  Blind  Test  Grid  S/N  map  is  shown 
in  Figure  6.  There  are  six  types  of  UXO  present  at  the  APG  Blind  Test  Grid  -  25  mm,  37  mm,  60 
mm,  81  mm,  105  mm  (IF=indirect  fire)  projectiles,  and  105  mm  HEAT  mortar  (DF=direct  fire). 
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To  East  (m) 


Figure  6,  Signal-to-noise  ratio  map  of  the  APG  Blind  Test  Grid.  Colors  are  plotted  on  a 
logarithmic  scale  for  better  visibility.  The  detection  threshold  of  S/N=4  is  equal  to  0.6  on 
the  logarithmic  scale. 


Principal  polarizabilities  from  the  soundings  that  have  been  selected  for  the  discrimination  were 
analyzed  using  a  variant  of  the  empirical  likelihood  ratio  method  (outlined  in  Gasperikova,  et  al. 
2009).  In  the  empirical  likelihood  ratio  method,  probability  densities  of  principal  polarizability 
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responses  were  estimated  from  the  training  data,  previously  eollected  UXO  and  scrap  responses 
separately. 

We  condensed  each  polarizability  response  down  to  a  set  of  nfeat  =10  numbers,  nine  of  which 
were  averages  of  products  of  major,  intermediate,  and  minor  principal  polarizabilities  with  time, 
in  three  time  bands.  The  logarithm  of  their  vector  magnitude  was  used  as  the  tenth  parameter, 
and  the  first  nine  values  were  normalized  by  that  magnitude.  Parameter  vectors  of  this  form  were 
differenced  with  the  vector  of  their  component  median  values  over  the  training  data  for  their 
class  (UXO  or  scrap)  for  use  in  forming  damped  trimmed  covariance  matrices  and  for 
UXO  and  scrap  responses,  which  were  used  in  forming  the  empirical  probability  densities  for 
UXO  and  scrap  training  data. 


The  empirical  densities  for  polarizability  responses  of  UXO  and  scrap  were  formed  by  centering 
a  generalized  Cauchy  density  at  each  of  the  different  polarizability  responses  from  the  class  and 
summing  over  the  responses  from  the  class.  The  generalized  Cauchy  density  that  was  used  for  a 
class  was  based  on  the  damped  trimmed  estimate  of  the  response  covariance  matrix  for  the  class 
Cciass^  with  length  scale  reduced  by  a  factor  of  (n‘^’^®®)'^^”feat,  so  that  the  contribution  from  each 
response  filled  of  the  volume  occupied  by  the  whole  distribution,  where  was  the 

number  of  responses  being  summed  over  to  form  the  density  estimate.  The  overall  form  was 


p(class)  (class) )  -  ^ 


iin  class 


(class)  , 


\  _|_  )*  ^y(class)  _  y(class)^  ^j^(class) 


where  y  =  0.2986/  ,  and  vj  was  the  vector  point  where  the  density  was  being  evaluated,  Vi 


were  the  training  response  vectors  the  density  was  based  on,  and  K  was 
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where  r()  is  the  gamma  function,  and  was  a  parameter  which  was  adjusted  using  cross 
validation  to  maximize  the  likelihood  of  the  training  data. 


The  probability  that  a  response  vector  vj  is  due  to  UXO  is 


(uxo) 


a„ 


aLf'"”’(vj)  +  aLpf'‘”'’'(Vj) 


2  2 

where  a  uxo  and  a  scrap  are  a  priori  relative  probabilities  of  UXO  and  scrap,  and  only  the  ratio  of 
a\xo/a^scrap  is  significant. 


In  artillery  ranges  where  no  cleanup  has  been  done,  assuming  a  dud  rate  of  0.1  and  that  ordnance  is 
blown  into  five  detectable  pieces,  a  reasonable  value  for  a  uxo/a  scrap  would  be  0.02.  However,  in  a 
small  test  grid,  constructed  to  emphasize  the  ability  of  equipment  to  discriminate  between  UXO 
and  scrap,  it  is  reasonable  to  emplace  a  much  higher  ratio  of  UXO  to  scrap,  such  as  1:1,  which 
would  suggest  using  a  value  of  a  uxo/a  scrap  =  1.  In  our  2006  survey  at  the  Yuma  Proving  Ground 
(YPG)  Blind  Test  Grid  we  estimated  there  were  230  occupied  cells  and  183  UXO,  yielding  a 
UXO/scrap  ratio  of  183/47  suggesting  using  a  uxo/a  scrap  =  3.89.  Using  this  value,  all  responses 
with  p^“™^(Vj)>  0.  0005  were  considered  UXO.  Probabilities  estimated  using  the  empirical 

likelihood  ratio  discrimination  method  tend  to  vary  by  orders  of  magnitude  so  the  exact  value  of 
a\xo/a^scrap  Only  affected  four  to  six  objects  with  near  the  cut-off  value. 
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Just  as  UXO  and  scrap  class  polarizability  response  densities  were  estimated,  we  estimated 


densities  for  the  polarizability  responses  of  different  types  of  UXO.  However,  due  to  the  limited 
number  of  responses  of  individual  UXO  in  our  response  library,  we  used  the  overall  UXO 
trimmed  covariance  matrix  estimate  in  the  place  of  individual  UXO  type  covariance 
matrices  in  their  estimated  distributions.  Similarly,  for  individual  UXO  types,  in  the  generalized 
Cauehy  density  formula  for  the  exponent  we  used  the  exponent  estimated  for  the  overall 
distribution  of  UXO  responses  m“’‘°.  Given  empirieal  distributions  for  elasses  elass  k, 

and  a  priori  expectations  of  relative  frequencies  a  ciass_k  for  the  different  classes,  the  probability 
that  a  response  is  from  an  object  in  class  q  is 

2  r  (class  _q)/  x 

(class _q)/  \  _  ^class  q^ _ 

k 


Here,  serap  was  included  as  one  of  the  classes  in  the  denominator.  For  want  of  prior  expeetations 
as  to  the  expected  relative  frequeney  of  different  UXO  types,  the  different  ^  were  taken  to 

be  equal  for  all  the  classes  except  for  serap,  and  the  sum  of  a  ciass_k  over  UXO  classes  was  equal 
2 

to  the  above. 
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5.  PERFORMANCE  ASSESSMENT 


For  the  hand-held  prototype  performanee  evaluation  we  submitted  two  tables  with  following 
information: 

(a)  eell  is  empty  or  oeeupied 

(b)  objeet  depth  below  the  system 

(e)  identification  as  a  scrap  or  intact  UXO 

(d)  if  cell  was  identified  as  containing  UXO,  type  of  UXO 

(e)  cells  with  data  that  can’t  be  analyzed. 

The  first  approach  was  a  statistical  approach  described  in  Chapter  4.  The  second  approach  was  a 
template  match  approach.  In  both  approaches  the  cell  was  considered  empty  if  S/N  was  less  than 
four,  and  “can’t  analyze”  category  if  S/N  was  between  four  and  six. 


5.1  Statistical  Approach  -  Empirical  Likelihood  ratio  Method 

The  statistical  approach  used  reprocessed  data,  and  it  is  described  in  Chapter  4.  In  this  approach 
from  400  fiag/cell  locations,  eight  were  labeled  as  “can’t  analyze”,  and  178  were  identified  as 
empty,  so  the  discrimination  was  done  on  214  flags.  From  those,  32  were  identified  as  scrap  and 
safe  to  leave  in  the  ground,  and  182  were  identified  as  UXO.  This  was  based  on  the  probability 
of  being  UXO  (p(UXO))  larger  or  equal  to  0.0005.  If  one  would  use  p(UXO)>0.0015,  the 
number  of  scrap  metals  would  increase  to  36.  Figure  7  shows  the  empty  cells  in  green,  occupied 
cells  in  red,  and  “can’t  analyze”  locations  in  blue.  Figure  8  shows  the  location  of  UXO  in  red  and 
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To  North  (m) 


scrap  in  green.  Figure  9  shows  the  recovered  object  depths.  All  objects  in  the  Blind  Test  Grid 


were  shallower  than  0.5  m. 


APG  Blind  Test  Grid  -  Statistical 


T 


^  empty 
^  occupied 
▲  can't  analyze 


A 


Figure  7,  The  APG  Blind  Test  Grid  cells  identified  as  empty  (green  circle),  occupied  (red 
diamond)  or  “can’t  analyze”  (blue  triangle)  using  the  statistical  approach. 
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To  North  (m) 


APG  Blind  Test  Grid  -  Statistical 


Figure  8,  The  APG  Blind  Test  Grid  cells  identified  as  scrap  (green  circle)  or  UXO  (red 
diamond)  using  the  statistical  approach. 
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APG  Blind  Test  Grid  -  Statistical 


Figure  9,  The  APG  Blind  Test  Grid  -  recovered  object  depths. 


The  scoring  results  are  shown  in  Table  2.  In  this  approach  we  correctly  identified  96%  and  94% 
of  munitions  in  the  response  and  discrimination  stages,  respectively.  The  objects  missed  were 
some  of  60  mm,  81  mm,  and  105  mm  projectiles/mortars  at  the  depth  greater  than  eight  times  of 
the  objects  diameter.  98%  of  scrap  was  correctly  identified  in  the  response  stage,  and  the 
probability  of  false  positive  in  the  discrimination  stage  was  0.77.  The  background  alarm  rates  for 
the  response  and  discrimination  stages  were  0.10  and  0.07,  respectively.  The  efficiency  rate  at 
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the  operating  point  was  0.99,  with  the  false  positive  rejeetion  rate  of  0.21  and  the  baekground 


alarm  rejeetion  rate  of  0.32.  For  the  ease  with  no  loss  in  the  UXO  deteetion  (the  effieieney  rate 
equal  to  1.0),  both  false  positive  and  baekground  alarm  rejeetion  rates  were  zero;  henee  all 
deteeted  objeets  would  have  to  be  dug  out.  Table  2b  shows  the  Blind  Test  Grid  eorreet  UXO  type 
elassifieation  of  targets  that  were  eorrectly  diseriminated  as  munitions.  Overall  67%  of  munitions 
were  elassified  eorreetly. 


Table  2a,  Blind  Grid  scoring  results  for  the  statistical  approach 


1  Response  Stage 

Discrimination  Stage 

Munitions^ 

Scores 

P/'“:  by  type 

All  Types 

105-mm 

81/60-mm 

37/25-mm 

All  Types 

105-mm 

81/60-mm 

37/25-mm 

0.98 

0.96 

1.00 

1.00 

0.97 

0.94 

1.00 

1.00 

0.96 

0.90 

0.97 

1.00 

0.94 

0.87 

0.97 

1.00 

0.91 

0.79 

0.88 

0.93 

0.90 

0.75 

0.88 

0.93 

1  By  Depth^  \ 

0  to  4D 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

4D  to  8D 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

8D  to  12D 

0.56 

0.50 

0.00 

1.00 

0.44 

0.33 

0.00 

1.00 

Clutter 

Scores 

Pcd 

PfP 

1  By  Mass  \ 

By  Depth^ 

All  Mass 

0  to  0.25  kg 

>0.25  to 

1  kg 

>1  to 

10  kg 

All  Mass 

0  to  0.25  kg 

>0.25  to 

1  kg 

>1  to 

10  kg 

All  Depth 

0.99 

0.82 

0.98 

0.95 

1.00 

1.00 

0.77 

0.71 

0.79 

1.00 

0.95 

0.71 

0  to  0.15  m 

0.98 

0.96 

1.00 

1.00 

0.77 

0.74 

0.78 

1.00 

0.15  to  03  m 

0.94 

0.80 

1.00 

1.00 

0.75 

0.40 

0.86 

1.00 

0.3  to  0.6  m 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

1  Background  Alarm  Rates  | 

1 

1 

“The  two  numbers  to  the  right  of  the  all  types  munitions  result  are  an  upper  and  lower  90  pereent 
eonfidenee  level  for  an  assumed  binomial  distribution. 

’’All  depths  are  measured  to  the  eenter  of  the  objeet. 
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Table  2b,  Blind  Grid  correct  UXO  type  classification  of  targets  correctly  discriminated  as 
munitions  for  tbe  statistical  approach 


Size 

Percentage 

25mm 

100% 

37mm 

100% 

60mm 

100% 

81mm 

13% 

105mm 

67% 

105  artillery 

20% 

Overall 

67% 

5.2  Template-match  Approach 

The  template  match  approach  used  data  processed  in  the  field  in  the  first  pass  (no  repeat  or  cross 
measurements),  and  compared  polarizability  curves  of  unknown  objects  to  the  polarizabilities  of 
the  objects  in  the  database/library.  Parameters  used  in  this  approach  were  (a)  a  similarity  of 
major  polarizability  curves,  (b)  the  similarity  of  medium  and  minor  polarizability  curves,  and  (d) 
the  distance  between  major  and  medium/minor  polarizabilities.  In  this  approach  from  400 
flag/cell  locations,  18  were  labeled  as  “can’t  analyze”.  Fourteen  of  these  responses  had  S/N 
between  four  and  six,  three  locations  with  S/N  less  than  four  but  with  inconsistent  responses,  and 
one  with  a  relatively  high  S/N  but  again  with  an  inconsistent  response.  181  cells  with  S/N  less 
than  four  were  considered  empty,  and  therefore  the  discrimination  was  done  on  201  flags.  From 
those  54  were  identified  as  scrap,  and  safe  to  leave  in  the  ground,  and  147  were  identified  as 
UXO  or  objects  that  have  to  be  dug. 
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APG  Blind  Test  Grid  -  Template  Match 
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Figure  10,  The  APG  Blind  Test  Grid  cells  identified  as  empty  (green  circle),  occupied  (red 
diamond)  or  “can’t  analyze”  (blue  triangle)  using  the  template  match  approach. 
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Figure  11,  The  APG  Blind  Test  Grid  cells  identified  as  scrap  (green  circle)  or  UXO  (red 
diamond)  using  the  template  match  approach. 


The  scoring  results  are  shown  in  Table  3.  In  this  approach  we  correctly  identified  97%  and  97% 
of  munitions  in  the  response  and  discrimination  stages,  respectively.  We  were  not  able  to  detect 
some  of  105  mm  projectiles/mortars  at  the  depth  greater  than  eight  times  of  the  objects  diameter. 
We  correctly  discriminated  all  detected  munitions.  We  also  correctly  identified  97%  of  scrap  in 
the  response  stage,  and  the  probability  of  false  positive  in  the  discrimination  stage  was  0.49.  The 
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background  alarm  rates  for  the  response  and  discrimination  stages  were  0.08  and  0.05, 
respectively.  The  efficieney  rate  at  the  operating  point  was  1.0,  with  the  false  positive  rejection 
rate  of  0.49  and  the  background  alarm  rejection  rate  of  0.44.  Table  3b  shows  the  Blind  Test  Grid 
correct  UXO  type  classification  of  targets  correctly  discriminated  as  munitions.  Overall  72%  of 
munitions  were  classified  correetly.  In  this  case,  we  did  not  separate  105  mm  projectiles  from 
105  mm  mortars. 


Table  2a,  Blind  Grid  scoring  results  for  the  template  match  approach 


1  Response  Stage 

Discrimination  Stage  | 

Munitions^ 

Scores 

All  Types 

105-mm 

81/60-mm 

37/25-mm 

All  Types 

105-mm 

81/60-mm 

37/25-mm 

0.99 

0.96 

1.00 

1.00 

0.99 

0.96 

1.00 

1.00 

0.97 

0.90 

1.00 

1.00 

0.97 

0.90 

1.00 

1.00 

0.93 

0.79 

0.93 

0.93 

0.93 

0.79 

0.93 

0.93 

1  By  Depth^  \ 

0  to  4D 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

4D  to  8D 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

8D  to  12D 

0.67 

0.50 

1.00 

1.00 

0.67 

0.50 

1.00 

1.00 

Clutter 

Scores 

Pea 

PfP 

1  By  Mass  \ 

By  Depth^ 

All  Mass 

0  to  0.25  kg 

■■ 

All  Mass 

0  to  0.25  kg 

■■ 

All  Depth 

0.99 

0.55 

0.97 

0.93 

1.00 

1.00 

0.49 

0.24 

0.67 

1.00 

0.93 

0.43 

0  to  0.15  m 

0.97 

0.94 

1.00 

1.00 

0.46 

0.23 

0.67 

1.00 

0.15  to  03  m 

0.94 

0.80 

1.00 

1.00 

0.69 

0.40 

0.71 

1.00 

0.3  to  0.6  m 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

1  Background  Alarm  Rates  | 

“The  two  numbers  to  the  right  of  the  all  types  munitions  result  are  an  upper  and  lower  90-pereent 
confidenee  interval  for  an  assumed  binomial  distribution. 

’’All  depths  are  measured  to  the  center  of  the  object. 
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Table  3b,  Blind  Grid  correct  UXO  type  classification  of  targets  correctly  discriminated  as 
munitions  for  tbe  template  match  approach 


Size 

Percentage 

25mm 

93% 

37mm 

100% 

60mm 

33% 

81mm 

53% 

105mm 

77% 

Overall 

72% 

6.  CONCLUSIONS 

This  survey  showed  that  the  same  diserimination  eapabilities  afforded  by  the  eart-mounted 
system  are  available  in  the  hand-held  unit,  although  with  a  slightly  redueed  depth  of  deteetion.  If 
only  single  objeets  are  present,  taking  the  measurements  0.15  m  before  and  after,  and  at  the  flag 
loeation,  along  the  survey  lines,  resulted  in  less  than  10%  of  loeations  that  eould  not  be  analyzed, 
while  eorreetly  diseriminating  all  munitions  with  the  S/N  above  the  threshold.  Although  only  7% 
of  the  soundings  used  in  the  diserimination  were  from  loeations  different  from  the  flag  loeations, 
having  three  measurements  per  flag  helped  with  the  data  quality  eontrol,  and  they  were  essential 
espeeially  when  a  person  other  than  the  one  interpreting  the  data  aequired  the  data.  This  is  a  very 
important  finding,  sinee  it  simplifies  the  field  data  aequisition  and  makes  surveys  more  effieient. 
If  flag  loeations  were  to  eontain  multiple  objeets,  however,  multiple  measurements  around  the 
flag  would  be  required  to  identify  them  eorreetly. 
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9.  ACRONYMS 


APG 

Aberdeen  Proving  Ground 

BUD 

Berkeley  UXO  Diseriminator 

ESTCP 

Environmental  Seeurity  Teehnology  Certifieation  Program 

FPGA 

Field  Programmable  Gate  Array 

LBNL 

Eawrenee  Berkeley  National  Eaboratory 

RFS 

Richmond  Field  Station 

SERDP 

Strategic  Environmental  Research  and  Development  Program 

S/N 

Signal-to-noise 

uxo 

Unexploded  Ordnance 

YPG 

Yuma  Proving  Ground 
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TABLE  1.  The  Calibration  Grid  discrimnation  results 


Celll  p(UXO) 

obj depth(m) 

LBLJD 

ArmylD 

true depth(m) 

filename 

A1  I  0.997886 

0.7155 

105mmHEAT 

155mmP 

0.7 

CAL  x1.85  y2  Run  0120.bin  diff.avg.lnv 

0.3241 

57mm 

57mmP 

0.32 

CAL  x6  y2  Run  0127.bin  diff.avg.lnv 

■WllBBBBill 

0.1559 

BLU 

BLU-26 

0.12 

CAL  x10.15  y2  Run  0132.bin  diff.avg.lnv 

■iMBBHlIcll 

0.2464 

105mmHEAT 

Clutter-Frag 

0.14 

CAL  x14  y2  Run  0136.bin  diff.avg.lnv 

I  E1  I  0.999997 

0.1758 

105mmlF 

105mm(IF) 

0.05 

CAL  x18  y2  Run  0140.bin  diff.avg.lnv 

lOKEB 

0.2123 

25mm 

25mm(DF) 

0.2 

CAL  x22  y2  Run  0144.bin  diff.avg.lnv 

I  G1  I  0.999995 

0.289 

37mm 

37mm(DF) 

0.17 

CAL  x25.85  y2  Run  0147.bin  diff.avg.lnv 

■StlllliliTigEBl 

0.3767 

2.75" 

60mm(IF) 

0.34 

CAL  x30.15  y2  Run  0153.bin  diff.avg.lnv 

■n^Hi 

0.4788 

2.15" 

81mm(IF) 

0.43 

CAL  x34.15  y2  Run  0157.bin  diff.avg.lnv 

■heeh 

0.4838 

57mm 

105mm(DF) 

0.58 

CAL  x38  y2  Run  0160.bin  diff.avg.lnv 

I  K1  I  0.999805 

0.5494 

57mm 

105mm(IF) 

0.55 

CAL  X41.85  y2  Run  0163.bin  diff.avg.lnv 

0.256 

2.15" 

2. 75  Rocket 

0.18 

CAL  x1.85  y6  Run  0005.bin  diff.avg.lnv 

I  B2  I  0.999977 

0.3574 

20mm 

20mmP 

0.23 

CAL  x6  y6.multl.llst.2dlp.lnv.2 

isaitiiiiTiTaaa 

0.199 

BLU 

BLU-26 

0.14 

CAL  x9.85  y6  Run  0013.bin  diff.avg.lnv 

I  D2  I  0.154639 

0.1602 

105mmHEAT 

Clutter-Frag 

0.02 

CAL  x14  y6  Run  0018.bin  diff.avg.lnv 

laEEEB 

0.1695 

81mm 

81mm(IF) 

0.11 

CAL  x18  y6  Run  0022.bin  diff.avg.lnv 

I  F2  I  0.999977 

0.259 

25mm 

25mm(DF) 

0.19 

CAL  x22.15  y6  Run  0027.bin  diff.avg.lnv 

laaitiiiiTigEEi 

0.2546 

37mm 

37mm(DF) 

0.24 

CAL  x26  y6  Run  0030.bin  diff.avg.lnv 

I  H2  I  0.999999 

0.2684 

60mm 

60mm(IF) 

0.33 

CAL  x30  y6  Run  0034.bin  diff.avg.lnv 

0.5801 

60mm 

81mm(IF) 

0.49 

CAL  x33.85  y6  Run  0037.bin  diff.avg.lnv 

0.5584 

57mm 

105mm(DF) 

0.75 

CAL  x38.15  y6  Run  0395.bin  diff.avg.lnv 

0.7871 

57mm 

105mm(IF) 

0.69 

CAL  x42  y6. multi. list. 2dlp.lnv. 2 

I  A3  I  0.999997 

0.112 

37mm 

40mmP 

0 

CAL  x1.85  ylO  Run  0081.bin  diff.avg.lnv 

■ilcMlIiliTigEBl 

0.3676 

BDU 

BDU-28 

0.23 

CAL  x6  ylO. multi. list. 2dip.llm2.lnv. 2 

laclBaiWgiEi 

0.0482 

25mm 

Clutter-Frag 

0.11 

CAL  xIO  ylO  Run  0090.bin  diff.avg.lnv 

■HclMKiWglEi 

0.1937 

M42 

Clutter-Frag 

0.15 

CAL  x14  ylO  Run  0094.bin  diff.avg.lnv 

I  E3  I  0.999999 

0.1281 

37mm 

37mm(DF) 

0.09 

CAL  x18.85  ylO  Run  0100.bin  diff.avg.lnv 

0.3257 

25mm 

25mm(DF) 

0.3 

CAL  x22  ylO  Run  0103.bin  diff.avg.lnv 

I  G3  I  0.998628 

0.5194 

37mm 

37mm(DF) 

0.34 

CAL  x26  y1 0. multi. list. 2dip.lim.inv. 2 

■5iciitiiiiyi33g 

0.5039 

60mm 

60mm(IF) 

0.55 

CAL  x30  ylO  Run  0112.bin  diff.avg.lnv 

I  I3  I  0.999965 

0.5361 

60mm 

81mm(IF) 

0.6 

CAL  x34  ylO  Run  0116.bin  diff.avg.lnv 

0.4319 

57mm 

105mm(DF) 

0.88 

CAL  x38  y1 0. multi. list. 2dip.inv. 2 

I  K3  I  0.999999 

0.5447 

105mmlF 

105mm(IF) 

0.5 

CAL  x42  ylO  Run  0124.bin  diff.avg.lnv 

0.1665 

37mm 

40mmP 

0.11 

CAL  x2.15  y14  Run  0158.bin  diff.avg.lnv 

0.1756 

BDU 

BDU-28 

0.17 

CAL  x6  y14  Run  0161.bin  diff.avg.lnv 

0.1845 

25mm 

Clutter-Frag 

0.2 

CAL  xIO  y14  Run  0165.bin  diff.avg.lnv 

0.2246 

25mm 

Clutter-Frag 

0.12 

CAL  x14  y14  Run  0169.bin  diff.avg.lnv 

0.0651 

25mm 

25mm(DF) 

0.04 

CAL  x18  y14  Run  0172.bin  diff.avg.lnv 

I  F4  I  0.999549 

0.3496 

25mm 

25mm(DF) 

0.3 

CAL  X21.85  y14  Run  0178.bin  diff.avg.lnv 

miEEEH 

0.3396 

37mm 

37mm(DF) 

0.39 

CAL  x26.15  y14  Run  0407.bin  diff.avg.lnv 

■  I  ll  1 1 1  M 

0.3666 

57mm 

60mm(IF) 

0.6 

CAL  x30.15  y14  Run  0188.bin  diff.avg.lnv 

Mmmmm 

0.6603 

105mmlF 

81mm(IF) 

0.62 

CAL  x34.15  y14  Run  0192.bin  diff.avg.lnv 

I  J4  I  0.999685 

0.5795 

105mmlF 

105mm(DF) 

0.61 

CAL  x37.85  y14  Run  0194.bin  diff.avg.lnv 

0.1079 

25mm 

105mm(IF) 

1.05 

CAL  x42.15  y14  Run  0441.bin  diff.avg.lnv 

I  A5  I  0.999655 

0.3825 

2.75" 

2.75Rocket 

0.34 

CAL  x2  y18  Run  0233.bin  diff.avg.lnv 

lililiiEEEEEE! 

0.2372 

20mm 

20mmP 

0.19 

CAL  x6  y18  Run  0237.bin  diff.avg.lnv 

I  C5  I  0.999939 

0.1913 

M42 

M42 

0.21 

CAL  xIO  y18  Run  0241.bin  diff.avg.lnv 

IHiliiWBgElcI 

0.1791 

60mm 

Clutter-Frag 

0.08 

CAL  x14  y18  Run  0245.bin  diff.avg.lnv 

I  E5  I  0.999999 

0.3152 

60mm 

60mm(IF) 

0.05 

CAL  x18.15  y18  Run  0253.bin  diff.avg.elev20.inv 

laKEB 

0.2405 

25mm 

25mm(DF) 

0.26 

CAL  x22  y18  Run  0256.bin  diff.avg.elev20.inv 

I  G5  I  0.999999 

0.2924 

37mm 

37mm(DF) 

0.3 

CAL  x26.15  y18  Run  0261.bin  diff.avg.elev20.inv 

■SlilltlililiEglcI 

0.2915 

60mm 

60mm(IF) 

0.28 

CAL  x30.15  y18  Run  0266.bin  diff.avg.lnv 

I  I5  I  0.999999 

0.4523 

60mm 

81mm(IF) 

0.53 

CAL  x34.15  y18  Run  0270.bin  diff.avg.lnv 

0.5495 

105mmlF 

105mm(DF) 

0.56 

CAL  x38  y18  Run  0273.bin  diff.avg.lnv 

0.5387 

105mmlF 

105mm(IF) 

0.48 

CAL  x42.15  y18  Run  0278.bin  diff.avg.lnv 

0.078 

57mm 

155mmP 

1.08 

CAL  x1.7  y22  Run  0422.bin  diff.avg.lnv 

I  B6  I  0.999987 

0.2824 

57mm 

57mmP 

0.32 

CAL  x6.15  y22  Run  0318.bin  diff.avg.lnv 

lasiitiiiiiiTngg 

0.2824 

M42 

M42 

0.23 

CAL  xIO  y22  Run  0321.bin  diff.avg.lnv 

I  D6  I  0.999905 

0.1737 

105mmlF 

Clutter-Frag 

0.04 

CAL  x14.15  y22  Run  0326.bin  diff.avg.lnv 

0.1666 

105mmHEAT 

105mm(DF) 

0.04 

CAL  x17.85  y22  Run  0328.bin  diff.avg.lnv 

I  F6  I  0.989512 

0.3119 

25mm 

25mm(DF) 

0.26 

CAL  x22  y22  Run  0333.bin  diff.avg.lnv 

0.2241 

37mm 

37mm(DF) 

0.18 

CAL  x26  y22  Run  0338.bin  diff.avg.lnv 

1  H6  1  0.999999 

0.3981 

60mm 

60mm(IF) 

0.34 

CAL  x30  y22  Run  0342.bin  diff.avg.lnv 

laEEHI 

0.379 

2.75" 

81mm(IF) 

0.35 

CAL  x34  y22  Run  0346.bin  diff.avg.lnv 

J6  0.999998 

0.5612 

105mmlF 

105mm(DF) 

0.47 

CAL  x37.85  y22  Run  0349.bin  diff.avg.lnv 

K6  0.999999 

0.3542 

105mmlF 

105mm(IF) 

0.32 

CAL  x42.15  y22  Run  0356.bin  diff.avg.lnv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X1 .85-y2-Run-01 20.bin-diff.avg.inv 
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APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


Time  (s) 


Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X1 0.1 5-y2-Run-01 32.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X1 4-y2-Run-01 36.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X1 8-y2-Run-01 40.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x22-y2-Run-0144.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x25.85-y2-Run-0147.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X30.1 5-y2-Run-01 53.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X34.1 5-y2-Run-01 57.bin-diff.avg.inv 


45 


Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x38-y2-Run-0160.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X41 .85-y2-Run-01 63.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x1.85-y6-Run-0005.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x6-y6. multi. Iist.2dip.inv.2 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x9.85-y6-Run-0013.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x14-y6-Run-0018.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x18-y6-Run-0022.bin-diff.avg.inv 
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APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x22.15-y6-Run-0027.bin-diff.avg.inv 

1 — I — I - 1 - 1 - 1 - 1 — I — I — I  I  I 

- maximum 

- medium 

minimum 


Position  (x  y  z):  (-0.1856  0.1021  0.259)  ^ 


Time  (s) 


Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x26-y6-Run-0030.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x30-y6-Run-0034.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x33.85-y6-Run-0037.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x38.15-y6-Run-0395.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x42-y6. multi. Iist.2dip.inv.2 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X1 .85-y1 0-Run-0081  .bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x6-y10. multi. Iist.2dip.lim2.inv.2 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X1 0-y1 0-Run-0090.bin-diff.avg.inv 


Time  (s) 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x14-y10-Run-0094.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X1 8.85-y1 0-Run-01  OO.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x22-y1 0-Run-01 03.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x26-y10. multi. Iist.2dip.lim.inv.2 
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APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


Time  (s) 


Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x34-y1 0-Run-01 1  G.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x38-y10. multi. Iist.2dip.inv.2 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x42-y1 0-Run-01 24.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X2.1 5-y14-Run-01 58.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x6-y1 4-Run-01 61  .bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X1 0-y1 4-Run-01 65.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X1 4-y1 4-Run-01 69.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X1 8-y1 4-Run-01 72.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X21 .85-y1 4-Run-01 78.bin-diff.avg.inv 


75 


Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X26.1 5-y1 4-Run-0407.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X30.1 5-y1 4-Run-01 88.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X34.1 5-y1 4-Run-01 92.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x37.85-y1 4-Run-01 94.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x42.15-y14-Run-0441.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x2-y18-Run-0233.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x6-y18-Run-0237.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X1 0-y1 8-Run-0241  .bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x14-y18-Run-0245.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X1 8.1 5-y1 8-Run-0253.bin-diff.avg.elev20.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x22-y18-Run-0256.bin-diff.avg.elev20.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X26.1 5-y1 8-Run-0261  .bin-diff.avg.elev20.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X30.1 5-y1 8-Run-0266.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X34.1 5-y1 8-Run-0270.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x38-y18-Run-0273.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X42.1 5-y1 8-Run-0278.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x1.7-y22-Run-0422.bin-diff.avg.inv 
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APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


Time  (s) 


Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X1 0-y22-Run-0321  .bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-X1 4.1 5-y22-Run-0326.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x17.85-y22-Run-0328.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x22-y22-Run-0333.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x26-y22-Run-0338.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x30-y22-Run-0342.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x34-y22-Run-0346.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x37.85-y22-Run-0349.bin-diff.avg.inv 
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Principal  Polarizabilities  (m  Is) 


APPENDIX  1.  The  Calibration  Grid  polarizability  responses 


CAL-x42.15-y22-Run-0356.bin-diff.avg.inv 
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